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Combination of Differential D*^ Cross-Section Measurements 
in Deep-Inelastic ep Scattering at HERA 


Abstract 


HI and ZEUS have published single-differential cross sections for inclusive D*^-meson 
production in deep-inelastic ep scattering at HERA from their respective final data sets. 
These cross sections are combined in the common visible phase-space region of photon 
virtuality > 5 GeV^, electron inelasticity 0.02 <y < 0.7 and the D*^ meson’s trans¬ 
verse momentum pt{D*) >1.5 GeV and pseudorapidity |77(D*)| < 1.5. The combination 
procedure takes into account all correlations, yielding significantly reduced experimental 
uncertainties. Double-differential cross sections d^a/dQ^dy are combined with earlier D*^ 
data, extending the kinematic range down to 2^ > 1.5 GeV^. Perturbative next-to-leading- 
order QCD predictions are compared to the results. 
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1 Introduction 


Measurements of open eharm produetion in deep-inelastic electrorQ-proton scattering (DIS) at 
HERA provide important input for stringent tests of the theory of strong interactions, quantum 
chromodynamics (QCD). Previous measurements lfTl-l20l have demonstrated that charm quarks 
are predominantly produced by the boson-gluon-fusion process, yg —)■ cc, whereby charm pro¬ 
duction becomes sensitive to the gluon distribution in the proton. Measurements have been 
obtained both from the HERA-I (1992-2000) and HERA-II (2003-2007) data-taking periods. 

Different techniques have been applied at HERA to measure open-charm production in DIS. 
The full reconstruction of D or mesons ifTT-l^ fTOl - n^fTSlfT^rT^I^ . the long lifetime of 
heavy flavoured hadrons [r7]-l9l [T2llT4l[T7llT9l or their semi-leptonic decays [IT3l are exploited. 
After extrapolation from the visible to the full phase space, most of these data have already 
been combined [l2n at the level of the reduced cross-sections and have provided a consistent 
determination of the charm contribution to the proton structure functions, a measurement of the 
charm-quark mass mc{mc) and improved predictions for W- and Z-production cross sections at 
the LHC. However, the extrapolation procedure requires theoretical assumptions, which lead to 
theoretical uncertainties comparable in size to the experimental uncertainties [1211 . Moreover, 
this combination was restricted to inclusive DIS variables only, such as the photon virtuality, 
2^, and the inelasticity, y. Alternatively, the measured cross sections can be combined directly 
in the visible phase space. In this case, dependences on the theoretical input are minimised 
and the charm production mechanism can be explored in terms of other variables. Such a 
combination, however, is possible only for data with the same final state, covering a common 
visible phase space. The recently published differential cross-section measurements by HI 
111511181 and ZEUS [|20l for inclusive -meson production fulfil this requirement. The analysis 
of fully reconstructed mesons also offers the best signal-to-background ratio and small 
statistical uncertainties. 

In this paper, visible -production cross sections Il6l[l5l[l8l|20l| at the centre-of-mass en¬ 
ergy i/s' = 318 GeV are combined such that one consistent HERA data set is obtained and 
compared directly to differential next-to-leading-order (NLO) QCD predictions. The combi¬ 
nation is based on the procedure described elsewhere ll2Tl - l24ll . accounting for all correlations 
in the uncertainties. This yields a significant reduction of the overall uncertainties of the mea¬ 
surements. The possibility to describe all measurements both in shape and normalisation with 
a single set of theory parameter values is also investigated and interpreted in terms of future 
theory improvements. 

The paper is organised as follows. In Section [2] the theoretical framework is briefly intro¬ 
duced that is used for applying phase-space corrections to the input data sets prior to combina¬ 
tion and for providing NLO QCD predictions to be compared to the data. The data samples used 
for the combination are detailed in Section[3]and the combination procedure is described in Sec¬ 
tion HI The combined single- and double-differential cross sections are presented in Section [5] 
together with a comparison of NLO QCD predictions to the data. 


2 Theoretical predictions 

The massive fixed-flavour-number scheme (LENS) [l25ll is used for theoretical predictions, since 
it is the only scheme for which fully differential NLO calculations ll26l are available. The cross- 
section predictions for D*^ production presented in this paper are obtained using the HVQDIS 

*In this paper, ‘electron’ is used to denote both electron and positron if not stated otherwise. 





program [12^ which provides NLO QCD (0(a^)) calculations in the 3-flavour FFNS for charm 
and beauty production in DIS. These predictions are used both for small phase-space corrections 
of the data, due to slightly different binning schemes and kinematic cuts, and for comparison to 
data. 

The following parameters are used in the calculations and are varied within certain limits to 
estimate the uncertainties associated with the predictions: 

• The renormalisation and factorisation scales are taken as jir = iif = +4m^. The 

scales are varied simultaneously up or down by a factor of two for the phase-space cor¬ 
rections where only the shape of the differential cross sections is relevant. For absolute 
predictions, the scales are changed independently to 0.5 and 2 times their nominal value. 

• The pole mass of the charm quark is set to me = 1.50 ±0.15 GeV. This variation also 
affects the values of the renormalisation and factorisation scales. 

• For the strong coupling constant the value ^{Mz) =0.105 ±0.002 is chosen 11211 
which corresponds to ^ (M 2 ) = 0.116 ± 0 . 002 . 

• The proton parton density functions (PDFs) are described by a series of FFNS variants 
of the HERAPDFl .0 set Il24ll at NLO determined within the HERAFitter If27l framework, 
similar to those used in the charm combination paper |[2T]| . Charm measurements were 
not included in the determination of these PDF sets. For all parameter settings used here, 
the corresponding PDF set is used. By default, the scales for the charm contribution to the 
inclusive data in the PDF determination were chosen to be consistent with the factorisa¬ 
tion scale used in HVQDIS, while the renormalisation scale in HVQDIS was decoupled 
from the scale used in the PDF extraction, except in the cases where the factorisation 
and renormalisation scales were varied simultaneously. As a cross check, the renormal¬ 
isation scales for both heavy- and light-quark contributions are varied simultaneously in 
HVQDIS and in the PDF determination, keeping the factorisation scales fixed. The re¬ 
sult lies well within the quoted uncertainties. The cross sections are also evaluated with 
3-flavour NLO versions of the ABM |[2^ and MSTW |f29l PDF sets. The differences 
are found to be negligible compared to those from varying other parameters, such that no 
attempt for coverage of all possible PDFs is made. 

The NLO calculation performed by the HVQDIS program yields differential cross sections 
for charm quarks. These predictions are converted to D*^-meson cross sections by applying 
the fragmentation model described in a previous publication ll^ . This model is based on the 
fragmentation function of Kartvelishvili et al. If30ll which provides a probability density func¬ 
tion for the fraction of the charm-quark momentum transferred to the meson. The function 

is controlled by a single fragmentation parameter, aK- Different values of ax [l2TI are used 
for different regions of the invariant mass, s, of the photon-parton centre-of-mass system. The 
boundary 51 = 70±40 GeV^ between the first two regions is one of the parameter variations. 
The boundary §2 = 324 GeV^ between the second and third region remains fixed. The model 
also implements a transverse-fragmentation component by assigning to the meson a trans¬ 
verse momentum, kr, with respect to the charm-quark direction ll^ . The following parameters 
are used in the calculations together with the corresponding variations for estimating the uncer¬ 
tainties of the NLO predictions related to fragmentation: 

• The fragmentation parameter Uk, the bin boundary si and the average kr are varied 
according to a prescription described elsewhere ETl . 



The fraction of charm quarks hadronising into D*'^ mesons is set to 

/(c ^ D"+) = /(c ^ D"-) = 0.2287 ±0.0056 [|3T]| . 


The small beauty eontribution to the signal needs a detailed treatment of the B hadron 
deeay to mesons and is therefore obtained from NLO QCD predietions for beauty hadrons 
eonvoluted with deeay tables to D*^ mesons and deeay kinematies obtained from EvtGen [[^ . 
The parameters for the ealeulations and the uneertainties are: 


• The renormalisation and factorisation scales iir = Bf = + 4m^ are varied in the 

same way as deseribed above for eharm. The variations are applied simultaneously for 
the ealeulation of the eharm and beauty eross-seetion uneertainties. 

• The pole mass of the heauty quark is set to = 4.75 ± 0.25 GeV. 

• The fragmentation model for heauty quarks is based on the Peterson et al. If3^ parametri- 
sation using = 0.0035 ± 0.0020 If^ . 

• The fraction of beauty hadrons decaying into D*^ mesons is set to 

=0.173 ±0.020 0511. 

• The proton structure is described by the same PDF set (3-flavour scheme) used for the 
charm cross-section predictions. 

The total theoretical uncertainties are obtained by adding all individual contributions in 
quadrature. 


3 Data samples for cross-section combinations 

The HI 06] - [38ll and ZEUS 09ll detectors were general purpose instruments which consisted 
of tracking systems surrounded by electromagnetic and hadronic calorimeters and muon detec¬ 
tors. The most important detector components for the measurements combined in this paper are 
the central tracking detectors (CTD) operated inside solenoidal magnetic fields of 1.16 T (HI) 
and 1.43 T (ZEUS) and the electromagnetic sections of the calorimeters. The CTD of HI If37l 
(ZEUS [|40ll ) measured charged particle trajectories in the polar angular rangeEl of 15° < 0 < 
165(164)°. In both detectors the CTDs were complemented with high-resolution silicon ver¬ 
tex detectors: a system of three silicon detectors for HI, consisting of the Backward Silicon 
Tracker [|4T|| . the Central Silicon Tracker fl43 and the Forward Silicon Tracker [|43ll . and the 
Micro Vertex Detector Il44l for ZEUS. For charged particles passing through all active layers 
of the silicon vertex detectors and CTDs, transverse-momentum resolutions of (y{pT)/PT — 
O.OOIpt/ ©0.015 (HI) and o{pt)/pt ^0.0029pr/©0.0081 ©0.0012/p7’ (ZEUS), with pj 
in units of GeV, have been achieved. 

Each of the central tracking detectors was enclosed by a set of calorimeters comprising an 
inner electromagnetic and an outer hadronic section. The HI calorimeter system consisted 

^In both experiments a right-handed coordinate system is employed with the Z axis pointing in the nominal 
proton-beam direction, referred to as “forward direction”, and the X axis pointing towards the centre of HERA. 
The origin of the coordinate system is dehned by the nominal interaction point in the case of HI and by the centre 
of the CTD in the case of ZEUS. 





Data set 

Kinematic range 

.if 

(pb-i) 

2" 

(GeV2) 


Pt{D*) 

(GeV) 

rim 

I 

HI HERA-II (medium Q^) 

m 

5: 100 

0.02 :0.70 

> 1.5 

-1.5 : 1.5 

348 

II 

HI D*^ HERA-II (high Q^) 

El 

100: 1000 

0.02 :0.70 

> 1.5 

-1.5 : 1.5 

351 

III 

ZEUS D*^ HERA-II 

m 

5 : 1000 

0.02 :0.70 

1.5 : 20.0 

-1.5: 1.5 

363 

IV 

ZEUS D*^ 98-00 

m 

1.5 : 1000 

0.02 :0.70 

1.5 : 15.0 

-1.5: 1.5 

82 


Table 1: Data sets used in the eombination. For eaeh data set the respeetive kinematie range 
and the integrated luminosity, , are given. 


of the Liquid Argon ealorimeter (LAr) [l45ll and the baekward lead-seintillator ealorimeter 
(SpaCal) [[38l while the ZEUS deteetor was equipped with a eompensating uranium-seintillator 
ealorimeter (CAL) [|46l . Most important for the analyses eombined in this paper is the eleetro- 
magnetie part of the ealorimeters whieh is used to identify and measure the seattered eleetron. 
Eleetromagnetic energy resolutions o{E)/E of O.ll/V^ (LAr) [|4Tll . 0.07/\/E (SpaCal) Il48l 
and 0.18/ y/E (CAL), with E in units of GeV, were aehieved. 

The Bethe-Heitler proeess, ep ^ ejp, is used by both experiments to determine the lumi¬ 
nosity. Photons originating from this reaction were detected by photon taggers at about 100 m 
downstream of the electron beam line. The integrated luminosities are known with a precision of 
3.2% for the HI measurements HISIJISH and of about 2% for the ZEUS measurements lf^20[|49H . 

Combinations are made for single- and double-differential cross sections. In Table [U the 
dataset^ used for these combinations are listed together with their visible phase-space regions 
and integrated luminosities. The datasets I-III are used to determine single-differential com¬ 
bined cross sections as a function of the meson’s transverse momentum, px{D*), pseu¬ 
dorapidity, ri(D*), and inelasticity, z(D*) = (E(D*} — pz(D*))/(2Eey), measured in the lab¬ 
oratory frame, and of and y. The variables E{D*), pz{D*) and Eg denote the energy 
of the D*^ meson, the Z component of the momentum of the meson and the incoming 
electron energy, respectively. Owing to beam-line modifications related to the HERA-II high- 
luminosity running [f50l the visible phase space of these cross sections at HERA-II is restricted 

to 2 ^ > 5 

GeV^, which prevents a combination with earlier cross-section measurements 
for which the phase space extends down to 2^ = 1.5 GeV^. 

In the case of the double-differential cross section, d^(7/dyd2^, the kinematic range can be 
extended to lower using HERA-I measurements [l4l|6l|T0l|. In order to minimise the use of 
correction factors derived from theoretical calculations, the binning scheme of such measure¬ 
ments has to be similar to that used for the HERA-II data. One of the HERA-I measurements, 
set IV of Table [B satisfies this requirement and is therefore included in the combination of 
this double-differential cross section. The visible phase spaces of the combined single- and 
double-differential cross sections are summarised in Tabled 

The measurements to be combined for the single- and double-differential cross sections are 
already corrected to the Bom level with a running fine-stmcture constant a and include both 
the charm and beauty contributions to production. The total expected beauty contribution 
is small, varying from ~ 1% at the lowest 2^ to ~ 7% at the highest Q^. The cross sections 


^Of the two sets of measurements in ifTSl . that compatible with the above cuts is chosen. 


















Range in 

single double 

differential eross seetion 

(GeV2) 

Pt{D^) (GeV) 
|77(D*)| 

5-1000 1.5-1000 

0.02-0.7 

> 1.5 

< 1.5 


Table 2: Visible phase spaee of the eombined eross seetions. 

measured previously [I6l[l5l[l8l| are here eorreeted to the PDG value [l^ of the D® branehing 
ratio. 


3.1 Treatment of data sets for single-differential cross sections 

In order to make the input data sets eompatible with the phase spaee quoted in Tableland with 
eaeh other, the following eorreetions are applied before the eombination: 


• The HI eollaboration has published measurements of eross seetions separately for 
5 GeV^ <Q^ < 100 GeV^ (set I) and for 100 GeV^ < < 1000 GeV^ (set II). Due 

to the limited statisties at high a eoarser binning in pt{D*), ri{D*), z{D*) and y was 
used in set II eompared to set I. Therefore the eross seetion in a bin i of a given observable 
integrated in the range 5 GeV^ < < 1000 GeV^ is ealeulated aeeording to 


a,(5 < eVGeV^ < 1000) = cj,(5 < Q^/GcV^ < 100) 


( 1 ) 


+ c7NL°(100 < < 1000) 


o-int(100 < eVGeV^ < 1000) 
c7NLO(ioo < Q^/G&N^ < 1000)' 


Here aint denotes the integrated visible eross seetion and stands for the NLO pre- 
dietions obtained from HVQDIS. In this ealeulation both the experimental uneertainties 
of the visible eross seetion at high and the theoretieal uneertainties as deseribed in 
Seetion [2] are ineluded. The eontribution from the region 100 GeV^ < < 1000 GeV2 

to the full 2^ range amounts to 4% on average and reaehes up to 50% at highest px{D*). 

• The bin boundaries used for the differential eross seetion as a funetion of differ be¬ 
tween sets I, II and set III. At low this is solved by eombining the eross-seetion 
measurements of the first two bins of set I into a single bin. For > 100 GeV2 no 
eonsistent binning seheme eould be defined direetly from the single-differential eross- 
seetion measurements. However, the measurements of the double-differential eross see- 
tions d^a/d2^dy have been performed in a eommon binning seheme. By integrating 
these eross seetions in y, single-differential eross seetions in are obtained also for 
Q^> 100 GeV^ whieh ean be used direetly in the eombination. 


• The eross-seetion measurements in set III are restrieted to pt{D*) < 20 GeV while there 
is no sueh limitation in the phase spaee of the eombination. Therefore these eross seetions 
are eorreeted for the eontribution from the range prip*) > 20 GeV using HVQDIS. This 
eorreetion is found to be less than 0.1%. 






3.2 Treatment of data sets for double-differential cross sections 

Since the restrietion to the same phase spaee in ^ does not apply for the eombination of the 
double-differential eross seetions in ^ and y, the HERA-I measurement, set IV, is also ineluded 
in the eombination. This allows an extension of the ki nematic range down to > 1.5 GeV^. 
The pt{D*) ranges of the measurements of sets III and IV are eorreeted in the same way as for 
the single-differential cross sections. 

To make the binning seheme of the HERA-I measurement compatible with that used for 
the HERA-II datasets, the binning for all datasets is revised. Cross seetions in the new bins 
are obtained from the original bins using the shape of the HVQDIS predictions as described in 
SeetionlU The new binning is given in Seetion[5](Table|9l). Bins are kept only if they satisfied 
both of the following eriteria: 

• The predicted fraetion of the eross seetion of the original bin eontained in the kinematieal 
overlap region in and y between the original and corrected bins is greater than 50% 
(in most bins it is greater than 90%). 

• The theoretical uncertainty from the correction procedure is obtained by evaluating all 
uneertainties diseussed in Section [2] and adding them in quadrature. The ratio of the 
theoretieal uneertainty to the uneorrelated experimental uneertainty is required to be less 
than 30%. 

This proeedure ensures that the effeet of the theoretieal uneertainties on the eombined data 
points is small. Most of the HERA-II bins are left unmodified; all of them satisfied the eriteria 
and are kept. Out of the 31 original HERA-I bins, 26 bins satisfy the eriteria and are kept. The 
data points removed from the eombination mainly eorrespond to the low-y region where larger 
bins were used for the HERA-I data. 


4 Combination method 


The eombination of the data sets uses the minimisation method developed for the eombina¬ 
tion of inelusive DIS eross seetions II22IJ24II . as implemented in the HERAverager program llSTll . 
Eor an individual dataset e the eontribution to the funetion is defined as 


Xsxp,ei^\^j) y~'. 


m 






■ m‘ 


\2 ■ 


( 2 ) 


Here /i'’® is the measured value of the eross seetion in bin i and 5,,e,stat and 5,,e,uncor are the 
relative eorrelated systematic, relative statistieal and relative uncorrelated systematic uncertain¬ 
ties, respectively, from the original measurements. The quantities m‘ express the values of the 
expected eombined eross seetion for eaeh bin i and the quantities bj express the shifts of the 
eorrelated systematie-uneertainty sourees j, in units of the standard deviation. Several data sets 
providing a number of measurements (index e) are represented by a total x^ funetion, whieh is 
built from the sum of the Xexp,e funetions of all data sets 

xLim\ bj) = Y,xLp,e{^\ bj) + Y,b ]. 


( 3 ) 




The combined cross sections m' are obtained by the minimisation of Xm with respect to m‘ and 

bj- 

The averaging procedure also provides the covariance matrix of the m' and the uncertain¬ 
ties of the y at the minimum. The bj at the minimum and their uncertainties are referred to 
as “shift” and “reduction”, respectively. The covariances V of the m‘ are given in the form 
V = Kncor + Lfc^sys The matrix Vuncor is diagonal. Its diagonal elements correspond to the 
covariances obtained in a weighted average performed in the absence of any correlated system¬ 
atic uncertainties. The covariance matrix contributions correspond to correlated systematic 
uncertainties on the averaged cross sections, such that the elements of a matrix V^y^ are obtained 

as (Vsys)o ^ given a vector of systematic uncertainties. It is worth noting that, 

in this representation of the covariance matrix, the number of correlated systematic sources is 
identical to the number of correlated systematic sources in the input data sets. 

In the present analysis, the correlated and uncorrelated systematic uncertainties are predom¬ 
inantly of multiplicative nature, i.e. they change proportionally to the central values. In equa¬ 
tion Q the multiplicative nature of these uncertainties is taken into account by multiplying the 
relative errors and 5,,e,uncor by the cross-section expectation mf In charm analyses the statis¬ 
tical uncertainty is mainly background dominated. Therefore it is treated as being independent 
of m'. For the minimisation of Xm iterative procedure is used as described elsewhere [|2^ . 

The 55 systematic uncertainties obtained from the original publications were examined for 
their correlations. Within each data set, most of the systematic uncertainties are found to be 
point-to-point correlated, and are thus treated as fully correlated in the combination. In to¬ 
tal there are 23 correlated experimental systematic sources and 5 theory-related uncertainty 
sources. A few are found to be uncorrelated and added in quadrature. For the combination 
of single-differential cross sections the uncorrelated uncertainties also include a theory-related 
uncertainty from the corrections discussed in Section[3l which varies between 0 and 10% of the 
total uncertainty and is added in quadrature. Asymmetric systematic uncertainties were sym¬ 
metrised to the larger deviation before performing the combination. Except for the branching- 
ratio uncertainty, which was treated as correlated, all experimental systematic uncertainties were 
treated as independent between the HI and ZEUS data sets. Since the distributions in pt{D*), 
r}{D*), z{D*), and y are not statistically independent, each distribution is combined sepa¬ 
rately. 


5 Combined cross sections 


The results of combining the HERA-II measurements ifTSlfT^llOl as a function of pt{D*), 
J] (D*), z(D*), and y are given in Tables [3]-|7l together with their uncorrelated and correlated 
uncertainties The total uncertainties are obtained by adding the uncorrelated and correlated 
uncertainties in quadrature. 

The individual data sets and the results of the combination are shown in Figures [U-[51 The 
consistency of the data sets as well as the reduction of the uncertainties are illustrated further by 
the insets at the bottom of Figures [H and jH The combinations in the different variables have a 
X^ probability varying between 15% and 87%, i.e. the data sets are consistent. The systematic 

^ A detailed breakdown of correlated uncertainties can be found on 

http://www.desy.de/hlzeus/dstar2015/ 





shift between the two input data sets is covered by the respective correlated uncertainties. The 
shifts and reductions of the correlated uncertainties are given in Table [8l The improvement of 
the total correlated uncertainty is due to small reductions of many sources. While the effective 
doubling of the statistics of the combined result reduces the uncorrelated uncertainties, the 
correlated uncertainties of the combined cross sections are reduced through cross-calibration 
effects between the two experiments. Typically, both effects contribute about equally to the 
reduction of the total uncertainty. 

The combined cross sections as a function of pt{D*), t] (D*), z(D*), and y are compared 
to NLO prediction^ in Figures - [TOl In general, the predictions describe the data well. 
The data reach an overall precision of about 5% over a large fraction of the measured phase 
space, while the typical theoretical uncertainty ranges from 30% at low to 10% at high Q^. 
The data points in the different distributions are statistically and systematically correlated. No 
attempt is made in this paper to quantify the correlations between bins taken from two different 
distributions. Thus quantitative comparisons of theory to data can only be made for individual 
distributions. 

In order to study the impact of the current theoretical uncertainties in more detail, the effect 
of some variations on the predictions is shown separately in Figure dH compared to the same 
data as in Figures |6l [8] and [TOl Only the variations with the largest impact on the respective 
distribution are shown in each case. 

1. The NLO prediction as a function of pt{D*) (Figure dH top) describes the data better by 
either 

• setting the charm-quark pole mass to 1.35 GeV or 

• reducing the renormalisation scale by a factor 2 or 

• increasing the factorisation scale by a factor 2. 

Simultaneous variation of both scales in the same direction would largely compensate and 
would therefore have a much smaller effect. 

2. The prediction for the z{D*) distribution (Figure dH bottom left) describes the shape 
of the data better if the fragmentation parameters are adjusted such that the boundary 
between the two lowest fragmentation regions ETl is varied from the default of 70 GeV^ 
to its lower boundary of 30 GeV^. 

3. The preference for a reduced renormalisation scale already observed for pt{D*) is con¬ 
firmed by the z{D*) distribution (Figure dH bottom right). However, the shape of the 
z{D*) distribution rather prefers variations of the charm mass and the factorisation scale 
in the opposite direction to those found for the pt{D*) distribution. The distributions of 
the other kinematic variables do not provide additional information to these findings [[53 . 

As stated before, within the large uncertainties indicated by the theory bands in Figures dl- 
[51 all distributions are reasonably well described. However, the above study shows that the dif¬ 
ferent contributions to these uncertainties do not only affect the normalisation, but also change 
the shape of different distributions in different ways. It is therefore not obvious that a variant of 

^The NLO QCD prediction for the beauty contribution to D** production, calculated as described in Section[2] 
can be found on 

http://www.desy.de/hlzeus/dstar2015/ 



the prediction that gives a good description of the distribution in one variable will also give a 
good description of the distribution in another. 

Based on the study if items 1.-3. above, a ‘customised’ calculation is performed to demon¬ 
strate the possibility of obtaining an improved description of the data in all variables at the same 
time, both in shape and normalisation, within the theoretical uncertainties quoted in Section [2l 
For this calculation, the following choices were made: 

• From the three options discussed in item 1. above, the second is chosen, i.e. the renor¬ 
malisation scale is reduced by a factor 2, with the factorisation scale unchanged. 

• The change of the fragmentation parameter = 30 GeV^, as discussed in item 2. above, 
is applied. 

• At this stage, the resulting distributions are still found to underestimate the data normal¬ 
isation. As the renormalisation and factorisation scales are recommended to differ by at 
most a factor of two [l53l . the only significant remaining handle arising from items 1. 
and 3. is the charm-quark pole mass. This mass is set to 1.4 GeV, a value which is also 
compatible with the partially overlapping data used for a previous dedicated study ll2TI of 
the charm-quark mass. 

• All other parameters, which have a much smaller impact If52l than those discussed above, 
are left at their central settings as described in Section [2l 

The result of this customised calculation is indicated as a dotted line in Figures - [TOl 
A reasonable agreement with data is achieved simultaneously in all variables. This a posteri¬ 
ori adjustment of theory parameters is not a prediction, but it can be taken as a hint in which 
direction theoretical and phenomenological developments could proceed. The strong improve¬ 
ment of the description of the data relative to the central prediction through the customisation 
of the renormalisation scale is in line with the expectation that higher-order calculations will 
be helpful to obtain a more stringent statement concerning the agreement of perturbative QCD 
predictions with the data. The improvement from the customisation of one of the fragmentation 
parameters and the still not fully satisfactory description of the z{D*) distribution indicate that 
further dedicated experimental and theoretical studies of the fragmentation treatment might be 
helpful. 

In general, the precise single-differential distributions resulting from this combination, in 
particular those as a function of pt{D*), ri{D*) and z{D*), are sensitive to theoretical and 
phenomenological parameters in a way which complements the sensitivity of more inclusive 
variables like and y. 

The combined double-differential cross sections with the uncorrelated, correlated and total 
uncertainties as a function of and y are given in Table HI The total uncertainty is obtained 
by adding the uncorrelated and correlated uncertainties in quadrature. The individual data sets 
as well as the results of the combination are shown in Figure [l2l Including data set IV slightly 
reduces the overall cross-section normalisation with respect to the combination of sets I-III 
only. The pull distribution of the combination is shown in Figure [l3l The combination has a 
probability of 84%, i.e. all data sets are consistent. The shifts and reductions of the correlated 
uncertainties are given in Table [8l 

detailed breakdown of correlated uncertainties can be found on 

http://www.desy.de/hlzeus/dstar2015/. 



These combined cross sections are compared to NLO prediction^ in Figure [HI The cus¬ 
tomised calculation is also shown. In general the predictions describe the data well. The data 
have a precision of about 5-10% over a large fraction of the measured phase space, while the 
estimated theoretical uncertainty ranges from 30% at low to 10% at high Q^. As well as the 
single-differential distributions, these double-differential distributions give extra input to test 
further theory improvements. 


6 Conclusions 


Measurements of -production cross sections in deep-inelastic ep scattering by the HI and 
ZEUS experiments are combined at the level of visible cross sections, accounting for their sys¬ 
tematic correlations. The data sets were found to be consistent and the combined data have 
significantly reduced uncertainties. In contrast to the earlier charm combination at the level of 
reduced cross sections, the present combination does not have significant theory-related uncer¬ 
tainties and in addition distributions of kinematic variables of the mesons are obtained. 
The combined data are compared to NLO QCD predictions. The predictions describe the data 
well within their uncertainties. Higher order calculations would be helpful to reduce the theory 
uncertainty to a level more comparable with the data precision. Further improvements in the 
treatment of heavy-quark fragmentation would also be desirable. 
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Pt{D*) 

da 

dpriD*) 

^uncor 

^cor 

^ot 

(GeV) 

(nb/GeV) 

(%) 

{%) 

(%) 

1.50 

1.88 

2.35 

6.4 

4.7 

8.0 

1.88 

2.28 

2.22 

4.9 

4.2 

6.4 

2.28 

2.68 

1.98 

3.7 

4.0 

5.5 

2.68 

3.08 

1.55 

3.5 

3.7 

5.1 

3.08 

3.50 

1.20 

3.7 

3.5 

5.1 

3.50 

4.00 

9.29 X 10-1 

3.2 

3.4 

4.7 

4.00 

4.75 

6.14 X 10-1 

3.0 

3.5 

4.6 

4.75 

6.00 

3.19 X 10-1 

3.1 

3.3 

4.5 

6.00 

8.00 

1.15 X 10-1 

3.8 

3.7 

5.3 

8.00: 

11.00 

3.32 X 10-2 

5.4 

3.7 

6.5 

11.00 

20.00 

3.80 X 10-3 

10.4 

6.4 

12.2 


Table 3: The eombined differential D*^-produetion eross seetion in the phase spaee given in 
Table l2l as a funetion of pt{D*), with its uneorrelated (5uncor), eorrelated (5cor) and total (5tot) 
uneertainties. 



da 

dri{D*) 

^ncor 

^cor 

^ot 



(nb) 

(%) 

{%) 

(%) 

-1.50 

-1.25 

1.36 

5.8 

4.3 

7.2 

-1.25 

-1.00 

1.52 

4.6 

4.0 

6.1 

-1.00 

-0.75 

1.59 

4.6 

4.0 

6.1 

-0.75 

-0.50 

1.79 

3.8 

3.5 

5.2 

-0.50 

-0.25 

1.83 

3.8 

3.3 

5.1 

-0.25 

0.00 

1.89 

3.8 

3.7 

5.3 

0.00 

0.25 

1.86 

4.0 

3.4 

5.2 

0.25 

0.50 

1.88 

4.0 

3.6 

5.4 

0.50 

0.75 

1.91 

4.1 

3.5 

5.4 

0.75 

1.00 

1.92 

4.3 

4.0 

5.9 

1.00 

1.25 

2.08 

4.7 

4.0 

6.1 

1.25 

1.50 

1.81 

6.3 

4.8 

7.9 


Table 4: The eombined differential -production cross section in the phase space given in 
Tabled as a function of t] (£>*), with its uneorrelated (5uncor), correlated (5cor) and total (5tot) 
uncertainties. 





































z(D-) 

da 

dz(£»*) 

(nb) 

^ncor 

(*) 

Ocor 

(%) 

^ot 

(%) 

0.00: 0.10 

3.28 

9.5 

5.9 

11.2 

0.10: 0.20 

7.35 

4.8 

6.3 

7.9 

0.20 : 0.32 

8.61 

3.5 

4.6 

5.7 

0.32 : 0.45 

8.92 

2.7 

3.9 

4.7 

0.45 : 0.57 

8.83 

1.8 

4.0 

4.3 

0.57 : 0.80 

4.78 

2.4 

5.1 

5.6 

0.80: 1.00 

0.63 

8.1 

10.2 

13.0 


Table 5: The eombined differential D*=*=-produetion eross seetion in the phase spaee given in 
Table |2] as a funetion of z{D*), with its uneorrelated (5uncor), correlated (5cor) and total (5tot) 
uncertainties. 


(GeV2) 

da 

dg^ 

(nb/GeV2) 

^ticor 

(%) 

Ocor 

(%) 

^o[ 

(%) 

5 : 8 

4.74 X 10^1 

4.0 

5.0 

6.4 

8 : 10 

2.96 X 10 1 

4.3 

3.8 

5.8 

10 : 13 

2.12 X 10-1 

3.8 

4.0 

5.6 

13 : 19 

1.24 X 10-1 

3.2 

3.8 

5.0 

19 : 28 

7.26 X 10-2 

3.5 

3.6 

5.0 

28 : 40 

3.97 X 10-2 

3.7 

4.0 

5.5 

40 : 60 

1.64 X 10-2 

4.4 

4.7 

6.4 

60 : 100 

7.45 X 10“2 

5.2 

3.9 

6.5 

100: 158 

2.08 X 10-2 

7.2 

5.3 

9.0 

158 : 251 

8.82 X 10-4 

7.6 

5.0 

9.1 

251 : 1000 

7.50 X 10-2 

12.0 

6.7 

13.3 


Table 6: The combined differential -production cross section in the phase space given in 
Table [2] as a function of Q^, with its uneorrelated (5uncor), correlated (5cor) and total (5tot) un¬ 
certainties. 































3^ 


da 

dv 

^ncor 

^cor 

^ot 



(nb) 

(%) 

(%) 

(%) 

0.02 : 

0.05 

12.1 

5.8 

9.1 

10.8 

0.05 : 

0.09 

18.8 

3.9 

4.6 

6.0 

0.09 : 

0.13 

17.0 

3.4 

4.3 

5.5 

0.13 : 

0.18 

13.4 

3.7 

4.2 

5.6 

0.18 : 

0.26 

11.2 

3.4 

3.7 

5.0 

0.26 : 

0.36 

7.65 

3.7 

4.2 

5.6 

0.36 : 

0.50 

4.78 

4.0 

5.3 

6.6 

0.50: 

0.70 

2.65 

5.6 

6.4 

8.5 


Table 7: The eombined differential D*^-produetion eross seetion in the phase spaee given in 
Table |2] as a funetion of y, with its uncorrelated (5uncor), correlated (5cor) and total (5tot) uncer¬ 
tainties. 















Data set 

Name 

I,II 

Hi CJC effieieney 

I,II 

Hi luminosity 

I,II 

Hi MC PDF 

I,II 

Hi eleetron energy 

I,II 

Hi eleetron polar angle 

I,II 

HI hadronie energy seale 

II 

HI fragmentation threshold at high 

I,II 

HI alternative MC model 

I,II 

HI alternative MC fragmentation 

I,II 

HI fragmentation threshold 

I 

HI high 2^ uneertainty 

III 

ZEUS hadronie energy seale 

III 

ZEUS eleetron energy seale 

III 

ZEUS Pt{^s) eorreetion 

III 

ZEUS M{Kn) window variation 

III 

ZEUS traeking effieieney 

III 

ZEUS b MC normafisation 

III 

ZEUS PHP MC normafisation 

III 

ZEUS diffraetive MC normafisation 

III 

ZEUS MC reweighting {pt{D*) and Q^) 

III 

ZEUS MC reweighting (ri(D*)) 

III 

ZEUS fuminosity (HERA-II) 

IV 

ZEUS fuminosity (98-00) 

I-IV 

Theory variation 

I-IV 

Theory Pf variation 

I-IV 

Theory variation 

I-IV 

Theory fongitudunaf frag, variation 

I-IV 

Theory transverse frag, variation 


Table 8: Sources of point-to-point correlated uncertainties. For each source the affected data 
sets are given, together with the shift (sh) and reduction factor (red) in the combination obtained 
after the first iteration. For sources which do not affect the combination of a given differential 
cross section, no shifts and reductions are quoted. 





























































































































































































































d^C7 

dQ^dy 

^ncor 

^cor 

5tOt 

(GeV2) 


(nb/GeV^) 

{%) 

(%) 

(%) 

1.5 : 3.5 

0.02 : 0.09 

4.76 

12.9 

2.5 

13.2 


0.09 : 0.16 

5.50 

11.3 

2.6 

11.5 


0.16 : 0.32 

3.00 

12.0 

2.6 

12.3 


0.32 : 0.70 

9.21 X 10 ^ 

20.5 

2.5 

20.7 

3.5 : 5.5 

0.02 : 0.09 

2.22 

11.3 

2.8 

11.6 


0.09 : 0.16 

1.98 

7.9 

2.7 

8.3 


0.16 : 0.32 

1.09 

20.2 

2.7 

20.4 


0.32 : 0.70 

3.47 X 10 ^ 

14.6 

2.6 

14.8 

5.5 : 9 

0.02 : 0.05 

1.06 

12.3 

4.4 

13.1 


0.05 : 0.09 

1.46 

7.8 

4.1 

8.8 


0.09 : 0.16 

1.32 

5.4 

4.3 

6.9 


0.16 : 0.32 

7.73 X 10 ^ 

4.9 

3.9 

6.3 


0.32 : 0.70 

2.51 X 10 ^ 

5.6 

4.2 

7.0 

9 : 14 

0.02 : 0.05 

5.20 X 10^^ 

13.0 

6.6 

14.6 


0.05 : 0.09 

7.68 X 10^1 

6.6 

3.9 

7.7 


0.09 : 0.16 

5.69 X 10 1 

4.6 

2.8 

5.4 


0.16 : 0.32 

4.12 X 10-1 

4.6 

3.1 

5.6 


0.32 : 0.70 

1.51 X 10-1 

5.6 

4.0 

6.9 

14: 23 

0.02 : 0.05 

2.29 X 10-1 

11.4 

6.3 

13.0 


0.05 : 0.09 

3.78 X 10-1 

6.5 

4.1 

7.7 


0.09 : 0.16 

2.90 X 10 1 

4.8 

3.3 

5.8 


0.16 : 0.32 

1.86 X 10 1 

5.0 

3.4 

6.0 


0.32 : 0.70 

6.92 X 10-2 

6.2 

4.4 

7.7 

23 : 45 

0.02 : 0.05 

6.91 X 10-2 

14.8 

8.2 

16.7 


0.05 : 0.09 

1.23 X 10-1 

5.9 

3.6 

6.9 


0.09 : 0.16 

1.14 X 10-1 

4.4 

3.0 

5.3 


0.16 : 0.32 

7.42 X 10-2 

4.3 

3.0 

5.2 


0.32 : 0.70 

3.21 X 10-2 

5.2 

3.7 

6.4 

45 : 100 

0.02 : 0.05 

6.16 X 10 ^ 

33.5 

11.1 

35.3 


0.05 : 0.09 

2.70 X 10-2 

11.0 

4.4 

11.8 


0.09 : 0.16 

2.05 X 10-2 

8.0 

3.7 

8.8 


0.16 : 0.32 

1.99 X 10-2 

5.4 

3.2 

6.3 


0.32 : 0.70 

7.84 X 10-3 

6.9 

4.0 

7.9 

100: 158 

0.02 : 0.32 

4.12 X 10-3 

8.2 

4.1 

9.2 


0.32 : 0.70 

2.18 X 10-3 

11.1 

4.1 

11.9 

158 : 251 

0.02 : 0.30 

1.79 X 10-3 

10.2 

4.4 

11.1 


0.30 : 0.70 

9.28 X 10-4 

11.6 

4.6 

12.5 

251 : 1000 

0.02 : 0.26 

1.31 X 10-4 

14.5 

4.7 

15.3 


0.26 : 0.70 

1.18 X 10-4 

12.7 

5.0 

13.6 


Table 9: The eombined double-differential D*^-produetion eross seetion in the phase spaee 
given in Table [2] as a funetion of and with its uneorrelated (5uncor), correlated (5cor) and 
total (5tot) uncertainties. 
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Figure 1: Differential -production cross section as a function of pt{D*). The open triangles 
and squares are the cross sections before combination, shown with a small horizontal offset for 
better visibility. The filled points are the combined cross sections. The inner error bars indicate 
the uncorrelated part of the uncertainties. The outer error bars represent the total uncertainties. 
The histogram indicates the binning used to calculate the cross sections. The bottom part shows 
the ratio of these cross sections with respect to the central value of the combined cross sections. 
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Figure 2: Differential -production cross section as a function of ri{D*). The open triangles 
and squares are the cross sections before combination, shown with a small horizontal offset for 
better visibility. The tilled points are the combined cross sections. The inner error bars indicate 
the uncorrelated part of the uncertainties. The outer error bars represent the total uncertainties. 
The histogram indicates the binning used to calculate the cross sections. 
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Figure 3: Differential D*^-produetion eross seetion as a function of z{D*). The open triangles 
and squares are the cross sections before combination, shown with a small horizontal offset for 
better visibility. The filled points are the combined cross sections. The inner error bars indicate 
the uncorrelated part of the uncertainties. The outer error bars represent the total uncertainties. 
The histogram indicates the binning used to calculate the cross sections. 
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Figure 4: Differential D*^-produetion eross seetion as a funetion of Q^. The open triangles 
and squares are the eross seetions before eombination, shown with a small horizontal offset for 
better visibility. The filled points are the eombined eross seetions. The inner error bars indieate 
the uneorrelated part of the uneertainties. The outer error bars represent the total uncertainties. 
The histogram indicates the binning used to calculate the cross sections. The bottom part shows 
the ratio of these cross sections with respect to the central value of the combined cross sections. 
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y 

Figure 5: Differential D*^-produetion eross seetion as a function of y. The open triangles and 
squares are the cross sections before combination, shown with a small horizontal offset for 
better visibility. The filled points are the combined cross sections. The inner error bars indicate 
the uncorrelated part of the uncertainties. The outer error bars represent the total uncertainties. 
The histogram indicates the binning used to calculate the cross sections. 
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Figure 6: Differential D*^-produetion eross section as a function of pt{D*). The data points 
are the combined cross sections. The inner error bars indicate the uncorrelated part of the 
uncertainties. The outer error bars represent the total uncertainties. Also shown are the NLO 
predictions from HVQDIS (including the beauty contribution) and their uncertainty band. A 
customised NLO calculation (dotted line, see text) is also shown. The bottom part shows the 
ratio of these cross sections with respect to the central value of the combined cross sections. 
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Figure 7: Differential -production cross section as a function of ri{D*). The data points 
are the combined cross sections. The inner error bars indicate the uncorrelated part of the 
uncertainties. The outer error bars represent the total uncertainties. Also shown are the NLO 
predictions from HVQDIS (including the beauty contribution) and their uncertainty band. A 
customised NLO calculation (dotted line, see text) is also shown. 
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Figure 8: Differential D*^-produetion eross seetion as a funetion of z{D*). The data points 
are the eombined eross seetions. The inner error bars indieate the uneorrelated part of the 
uneertainties. The outer error bars represent the total uneertainties. Also shown are the NLO 
predietions from HVQDIS (ineluding the beauty eontribution) and their uneertainty band. A 
eustomised NLO ealeulation (dotted line, see text) is also shown. 
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Figure 9: Differential -production cross section as a function of Q^. The data points are the 
combined cross sections. The inner error bars indicate the uncorrelated part of the uncertainties. 
The outer error bars represent the total uncertainties. Also shown are the NLO predictions from 
HVQDIS (including the beauty contribution) and their uncertainty band. A customised NLO 
calculation (dotted line, see text) is also shown. The bottom part shows the ratio of these cross 
sections with respect to the central value of the combined cross sections. 
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Figure 10: Differential D*^-produetion eross seetion as a funetion of y. The data points are the 
eombined eross sections. The inner error bars indicate the uncorrelated part of the uncertainties. 
The outer error bars represent the total uncertainties. Also shown are the NLO predictions from 
HVQDIS (including the beauty contribution) and their uncertainty band. A customised NLO 
calculation (dotted line, see text) is also shown. 
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Figure 11: Differential D*^-produetion eross seetion (ratio to data) as a funetion of pt{D*) 
(top) and z{D*) (bottom) eompared to NLO predietions with different variations: eharm-quark 
mass nic, renormalisation seale Pr, faetorisation seale Pf and fragmentation bin boundary 5i. 
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Figure 12: Double-differential D*^-produetion eross seetions as a funetion of ^ and y. The 
open triangles, squares and eireles are the eross seetions before eombination, shown with a 
small horizontal offset for better visibility. The filled points are the eombined eross seetions. 
The inner error bars indicate the uncorrelated part of the uncertainties. The outer error bars 
represent the total uncertainties. 
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Figure 13: The pull distribution for the eombination of the double-differential eross see- 
tions. 
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Figure 14: Double-differential D*^-produetion eross seetion as a funetion of and y. The 
data points are the combined cross sections. The inner error bars indicate the uncorrelated part 
of the uncertainties. The outer error bars represent the total uncertainties. Also shown are the 
NLO predictions from HVQDIS (including the beauty contribution) and their uncertainty band. 
A customised NLO calculation (dotted line, see text) is also shown. 



































































